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Constitutive  Model  for  Anisotropic  Creep  Behaviors 
of  Single-Crystal  Ni-Base  Superalloys  in  the  Low-Temperature, 
High-Stress  Regime 


YOON  SUK  CHOI,  TRIPLICANE  A.  PARTHASARATHY, 

CHRISTOPHER  WOODWARD,  DENNIS  M.  DIMIDUK,  and  MICHAEE  D.  UCHIC 

A  crystallographic  constitutive  model  is  developed  to  capture  orientation-sensitive  primary  and 
secondary  creep  behaviors  within  approximately  20  deg  from  the  [0  0  1]  orientation  in  single¬ 
crystal  superalloys  for  the  low-temperature  and  high-stress  regime.  The  crystal  plasticity-based 
constitutive  formulations  phenomenologically  incorporate  experimentally  observed  dislocation 
micromechanisms.  Specihcally,  the  model  numerically  delineates  the  nucleation,  propagation, 
and  hardening  of  a(112)  dislocations  that  shear  multiple  y'  precipitates  by  creating  extended 
stacking  faults.  Detailed  numerical  descriptions  involve  slip-system  kinematics  from  a/2(110) 
dislocations  shearing  the  y-phase  matrix,  a(112)  stacking  fault  dislocation  ribbons  shearing  the 
y'-phase  precipitate,  interactions  between  a/2(110)  dislocations  to  nucleate  a(112)  dislocations, 
and  interactions  between  the  two  types  of  dislocations.  The  new  constitutive  model  was 
implemented  in  the  finite-element  method  (FEM)  framework  and  used  to  predict  primary  and 
secondary  creep  of  a  single-crystal  superalloy  CMSX-4  in  three  selected  orientations  near  the  [0 
0  1]  at  1023  K  (750  °C)  and  750  MPa.  Simulation  results  showed  a  reasonable,  qualitative 
agreement  with  the  experimental  data.  The  simulation  results  also  indicated  that  a/2(110) 
matrix  dislocations  are  important  to  limit  the  propagation  of  a(112)  dislocations,  which  leads  to 
the  transition  to  secondary  creep. 
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I.  INTRODUCTION 

SlNGEE-CrySta.1  Ni-base  superalloys  have  been  an 
unbeatable  choice  for  hot-section  turbine  blades  because 
of  their  superior  thermomechanical  performance  at 
increased  temperatures  compared  with  other  high-tem- 
perature  materials.  It  is  generally  understood  that  a 
large  volume  fraction  of  coherent  cuboidal  Ni3Al-base 
(LI2)  y'-phase  precipitates  regularly  distributed  in  a 
Ni-base  solid  solution  y-phase  matrix  provides  an 
excellent  barrier  against  the  motion  of  dislocations  on 
thermomechanical  loading.  Experimental  studies^* 
show  that  the  dislocation  behavior  in  the  y  matrix  and 
the  y/y'  interface,  and  the  stability  of  y'  precipitates  vary 
depending  on  the  temperature  and  stress,  and  they  result 
in  different  creep  behaviors  for  these  alloys. 

Near-[0  0  1]  orientation  (<20  deg)  creep  behavior  of 
Ni-base,  single-crystal  superalloys  has  been  subject  to 
numerous  experimental  and  numerical  studies.  Experi- 
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mental  studies^^  identihed  three  distinctive  near-[0  0  1] 
creep  behaviors  depending  on  ranges  of  temperature  and 
stress.  Figure  1  illustrates  schematically  these  creep 
behaviors  along  with  corresponding  temperature-stress 
regimes  based  on  the  creep  data  for  a  second-generation 
Ni-base  single  crystal  superalloy  CMSX-4. Other 
second-generation,  single-crystal  superalloys,  such  as 
Rene  N5  and  PWA1484,  exhibit  similar  creep  behav- 
iors.^^^  In  Figure  1,  the  data  points  indicate  stresses  and 
temperatures  at  which  actual  creep  tests  were  taken. 
Here,  the  data  points  with  a  same  symbol  indicate  a 
similar  creep  behavior  (or  creep  curves)  as  grouped  by 
shaded  regions.  In  the  low-temperature,  high-stress 
regime  (EH  regime  in  Figure  1),  the  deformation  is 
dominated  by  primary  and  secondary  creep.  The  creep 
strain  and  time  for  the  transition  from  primary  to 
secondary  creep  is  highly  sensitive  to  the  crystallo¬ 
graphic  orientation.  The  magnitude  of  primary  creep  is 
often  tied  qualitatively  to  the  degree  of  single  slip  in  a 
particular  slip  system  until  it  is  hampered  by  latent  hard¬ 
ening  caused  by  activated  additional  slip  systems. 
Typically,  the  deformation-induced  geometrical  rota¬ 
tions  are  sufficient  to  activate  additional  slip  systems. 
In  the  intermediate-temperature,  intermediate-stress 
regime  (MM  regime  in  Figure  1),  a  typical  creep  curve 
shows  an  extended  period  of  small  creep  strain  (in  the 
absence  of  primary  creep)  until  it  reaches  a  tertiary- 
creep-like  stage. In  this  regime,  the  gradual  pro¬ 
gression  of  directional  coarsening  (e.g.,  rafting)  of 
y'  precipitates  is  believed  to  be  responsible  for  the 
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Fig.  1  Three  temperature  stress  regimes  grouped  in  a  stress  tern 
perature  plot:  Each  regime  represents  a  distinctive  creep  behavior,  as 
illustrated  schematically  a  creep  curve.  The  classification  of  creep 
behaviors  is  based  on  the  literature  survey  of  near  [0  0  1]  creep  of 
CMSX  4.  All  units  in  the  schematic  creep  curves  are  arbitrary. 


tertiary-dominated  creep  behavior  as  illustrated  in 
Figure  1.^*'  Creep  in  the  high-temperature,  low-stress 
regime  (HL  regime  in  Figure  1)  typically  exhibits  a  small 
amount  of  primary  creep  followed  by  clear  secondary 
creep  and  catastrophic  tertiary  creep  that  immediately 
causes  failure.  In  this  regime,  the  completion  of  y' 
directional  coarsening  is  observed  even  in  the  primary 
creep  stage,  and  the  progression  of  secondary  and 
tertiary  creep  is  believed  to  be  controlled  by  the 
destabilization  of  the  y/y'  interfacial  dislocation  network 
and  the  heterogeneity  of  inherent  defects. Note 
that  these  three  distinctive  creep  regimes  reflect  different 
microstructural  sensitivities  and  dislocation  micromech¬ 
anisms.  Thus,  it  is  imperative  to  incorporate  those 
micromechanisms  into  a  constitutive  model  to  capture 
microstructure  sensitivities  of  the  creep  behavior  com¬ 
prehensively  in  all  stress  and  temperature  regimes.  It  is 
the  authors’  long-term  goal  to  develop  a  comprehensive 
creep  constitutive  model  for  all  stress  and  temperature 
regimes. 

As  a  first  step,  in  the  current  study,  we  focused  on 
developing  a  constitutive  model  that  represents  the 
orientation-sensitive  primary  and  secondary  near-[0  0  1] 
creep  behavior  in  the  LH  regime.  Significant  progress  has 
been  achieved  in  understanding  dislocation  micromech¬ 
anisms  in  this  creep  regime, although  a  mechanistic 
strain-hardening  law  for  creep  has  not  been  achieved. 
Creep  models  have  also  improved  progressively  by  these 
pacing  advances  in  understanding  of  creep  mechanisms. 
The  first  part  of  this  paper  summarizes  the  current 
understanding  of  dislocation  micromechanisms  and  creep 
models  for  orientation-sensitive  creep  in  the  LH  regime. 
Distinctive  features  of  the  new  creep  model  developed 
in  the  current  study  are  discussed  subsequently  along 
with  detailed  constitutive  formulations.  In  particular. 


representative  microstructural  features  that  the  new  creep 
constitutive  model  accounts  for  are  clarified.  Finally, 
creep  at  1023  K  (750  °C)  and  750  MPa  was  simulated  for 
selected  loading  directions  (within  1 8  deg  from  the  [001]), 
and  the  results  are  compared  with  the  experimental  creep 
curves  for  CMSX-4. 


II.  DISLOCATION  MICROMECHANISMS 

AND  PREVIOUS  MODELS  FOR  NEAR-[0  0  1] 
CREEP  IN  THE  LH  REGIME 

In  the  LH  regime,  the  stress  is  high  enough  to  induce 
y'  shearing  while  the  temperature  is  still  too  low  to 
facilitate  directional  coarsening  of  /  precipitates.  Exper¬ 
imental  observations  indicated  that  dislocations  with  a 
net  Burgers  vector  a(112)  cut  y'  precipitates  by  creating 
stacking  faults  (SF  shearing)  in  a  form  of  a  superlattice 
intrinsic  stacking  fault  (SISF)  or  a  superlattice  extrinsic 
stacking  fault  (SESF)  in  this  temperature  and  stress 
regime.^  Here,  a/2(112),  a/3(112),  and  a/6(112) 
partial  dislocations  with  a  net  Burgers  vector  a(112) 
are  formed  in  various  combinations  to  shear  /  precip¬ 
itates. However,  current  constitutive  models  have  not 
been  developed  to  a  sufiicient  state  as  to  represent 
accurately  the  spectrum  of  possible  a(  112) -type  disloca¬ 
tion  reactions  and  the  full  fidelity  of  their  effects  on 
strain  hardening  during  creep,  nor  has  that  challenge 
been  solved  by  the  present  model.  In  the  current  work, 
these  partial  dislocation  reactions  are  referred  to  as 
“a(112)  dislocation  ribbons”  without  the  subsequent 
distinction  between  specific  dislocation  reactions.  Rae 
et  proposed  micromechanisms  for  primary  creep 

caused  by  SF  shearing  of  y'  precipitates  by  a(112) 
dislocation  ribbons  and  the  transition  to  secondary 
creep.  Their  primary  creep  mechanism  involves  three 
steps:  the  nucleation,  propagation,  and  hardening  of 
a(112)  dislocation  ribbons.  The  nucleation  of  a(112) 
dislocation  ribbons  requires  a  reaction  of  two  active  a/ 
2(110)  dislocations  in  the  y-matrix  channels  or  y/y' 
interfaces.  Once  a(112)  dislocation  ribbons  are  nucle¬ 
ated,  their  resolved  shear  stresses  should  be  high  enough 
to  percolate  dislocations  through  y'  precipitates  (in  a 
form  of  mixed  cutting  and  looping).  Shearing  of  y' 
precipitates  can  be  done  by  a  pair  of  partial  dislocations 
(of  a  net  Burgers  vector  a(112))  separated  by  SISF  or 
SESF,  or  two  pairs  of  partial  dislocations  (still  of  a  net 
Burgers  vector  a(112))  separated  by  the  complex  stack¬ 
ing  fault  and  the  antiphase  boundary  (APB).^^’^’’^ 
Generally,  the  a/2(110)  dislocations  are  not  observed 
to  cut  /  precipitates  because  an  APB  fault  would  have 
to  be  formed  and  then  annihilated  by  a  second  coupled 
a/2(110)  dislocation.  Here,  highly  orientation-sensitive 
primary  creep  in  the  vicinity  of  the  [0  0  1]  orientation  is 
believed  to  be  caused  by  the  propagation  of  single 
a(112)  {1  1  1}  slip.  Suppose  that  orientations  A  and  B 
are  located  in  slightly  different  positions  in  the  stereo¬ 
graphic  projection  near  the  [0  0  1],  as  shown  as  open 
symbols  in  Figure  2,  and  a(112){l  1  1}  slip  systems  are 
operative.  After  applying  tension  stress  in  A  and  B 
directions,  different  amounts  of  [112](111)  slip  (hence, 
primary  creep)  are  expected  for  orientations  A  and  B 
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Fig.  2  Stereographic  projection  near  [0  0  1],  The  arrows  indicate 
rotation  paths  of  orientations  A  and  B  once  the  [ll2](Ill)  slip  sys 
tern  is  activated  under  the  [0  0  1]  tension. 

because  the  duration  of  [112](111)  slip  may  be  different 
between  the  two  orientations  as  differing  amounts  of 
geometrical  rotation  (indicated  by  dotted  arrows  in 
Figure  2)  are  required  to  reach  a  multislip  orientation 
(orientations  along  the  [0  0  l]-[0  1  1]  symmetry  bound¬ 
ary  in  Figure  2).  A  simple  dislocation  percolation  based 
calcnlation  suggested  that  SF  shearing  of  y'  precipitates 
by  (112){1  1  1}  slip  produces  much  more  shear  strain 
(by  a  factor  of  approximately  13)  compared  with  the 
case  for  only  y-matrix  shearing  by  a/2(l  10){1  1  1}  slip.^"*^ 
That  estimate  indicates  that  a  relatively  small  difference 
in  single  a(112){l  1  1}  slip  duration  could  result  in  a 
large  difference  in  primary  creep  strain,  ffere,  single 
a(112){l  1  1}  slip  is  hardened  by  interactions  between 
a(112)  dislocation  ribbons  and  a/2(110)  matrix  disloca¬ 
tions,  as  well  as  between  a(112)  dislocation  ribbons 
from  different  a(112){l  1  1}  slip  systems. This  leads  to 
the  transition  to  secondary  creep. 

Now,  it  seems  to  be  clear  that  the  orientation-sensitive 
primary  creep  for  crystals  oriented  near  [0  0  1]  neces¬ 
sitates  relatively  easy  propagation  of  a(112)  dislocation 
ribbons  throughout  the  entire  y/y'  microstructure.  In 
other  words,  the  orientation  sensitivity  may  be  less 
pronounced  if  the  propagation  of  a(112)  dislocation 
ribbons  is  hampered  effectively.  This  is  probably  the 
reason  for  different  levels  of  orientation  sensitivity  of 
near-[0  0  1]  primary  creep  observed  for  a  variety  of 
single-crystal  superalloys. Different  compositions 
and  heat  treatments  result  in  different  values  of  fault 
energies  (such  as  SISF  and  APB  energies),  the  lattice 
misfit  between  y  and  y'  phases,  and  different  volume 
fractions  of  /  precipitates.  All  these  factors  influence  the 
efficiency  of  a(112)  dislocation  ribbon  propagation, 
hence  the  orientation  sensitivity  of  near-[0  0  1]  primary 
creep.  Unfortunately,  state-of-the-art  constitutive  mod¬ 
eling  approaches  cannot  capture  fully  such  chemistry 
and  microstructure  dependences.  The  current  study  is 
limited  to  representing  primary  and  secondary  creep 
behaviors  phenomenologically  for  selected  single-crystal 
superalloys  showing  highly  orientation-sensitive  near-[0 
0  1]  creep. 


Early  modeling  of  the  anisotropic  creep  behavior  near 
the  [0  0  1]  was  done  by  Matan  et  alP^  They  proposed  a 
phenomenological  creep  model,  which  was  adopted 
from  Gilman’s  dislocation  density  model, by  assum¬ 
ing  that  the  mobile  dislocation  density  directly  ties  to  the 
creep  strain.  To  capture  the  orientation  dependence  of 
near-[0  0  1]  creep,  they  set  the  mobile  dislocation 
attrition  coefficient  to  be  inversely  proportional  to  9  (the 
rotation  required  to  reach  a  multislip  orientation  from 
the  original  loading  orientation).  MacLachlan  et  alP^ 
proposed  a  series  of  creep  models  for  anisotropic  creep 
of  single-crystal  superalloys.  Their  first  model  was  a 
continuum  damage  model,  which  incorporated  slip 
system  kinematics  to  predict  the  anisotropic  creep 
behavior. They  introduced  plastic  shear  rates  and 
resolved  shear  stresses  for  a/2(011)  slip  and  a(112)  slip, 
and  their  effective  values  were  linked  to  continuum 
damage  parameters.  In  addition,  they  used  a  kinematic 
hardening  term  to  account  for  hardening  contributions 
caused  by  interactions  between  a/2(011)  and  a(112)  slip. 
Their  first  model  was  refined  even  more  by  enhancing 
softening  and  hardening  terms  for  a(112)  slip  and,  the 
modified  model  captured  the  near-[0  0  1]  anisotropic 
creep  behavior  of  CMSX-4  at  1023  K  (750  °C)  and 
750  MPa.^*®'  In  their  most  recent  model,  hardening 
contributions  from  a/2 (01 1)  slip  interactions,  a(l  12)  slip 
interactions,  and  a/2(011)-a(112)  slip  interactions  were 
treated  separately. ^  This  model  was  used  to  predict  the 
creep  behavior  of  a  single-crystal  turbine  blade  compo¬ 
nent. Recently,  Ma  et  proposed  an  analytical 
creep  model  for  the  prediction  of  anisotropic  creep  for 
CMSX-4.  They  treated  creep  strain  rates  in  y  matrices 
and  /  precipitates  by  a/2(011)  slip  and  a(112)  slip, 
respectively.  Their  creep  model  was  based  on  the 
evolution  of  dislocation  densities  in  y  matrices  and  / 
precipitates,  as  well  as  corresponding  dislocation  veloc¬ 
ities.  Major  microstructnral  parameters,  such  as  a  width 
of  y-matrix  channels  (a  volume  fraction  of  y'  precipi¬ 
tates),  the  APB  energy,  and  the  y/y'  lattice  misfit,  were 
incorporated  in  formulations  of  dislocation  densities 
and  velocities  but  not  tied  explicitly  to  slip  system 
activation.  They  also  incorporated  an  incubation  effect, 
which  is  often  observed  in  low-temperatnre,  relatively 
low-stress  creep,  and  the  effect  of  void  damage  for  the 
better  prediction  of  primary  and  tertiary  creep  behav¬ 
iors,  respectively. 

Generally,  constitutive  models  for  near-[0  0  1]  creep  in 
the  LH  regime  focused  on  the  facilitation  of  single 
a(112)  slip  in  the  primary  creep  stage  and  latent 
hardening  of  a(112)  slip  (a  result  of  the  single  slip- 
induced  geometrical  rotation  and  the  interaction  with 
a/2(110)  y-matrix  dislocations)  for  the  transition  from 
primary  to  secondary  creep.  In  crystal  plasticity  (CP) 
FEM,  microstructnral  influences  on  single  slip  and  strain 
hardening  in  a  single  crystal  may  be  less  well  represented 
because  the  single-slip  behavior  is  treated  in  a  continuum 
basis,  which  means  no  discrete  representation  of  dislo¬ 
cation  transmission  or  interaction  in  and  across  precip¬ 
itates  in  a  single  crystal.  Because  of  this,  constitutive 
models  often  used  “flags”  (in  a  form  of  hyperbolic  or 
sinusoidal  functions)  to  facilitate  and  control  the  single 
slip  behavior  in  CP-EEM  effectively.^'®  The  current 
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constitutive  model  also  used  those  functional  forms  for 
the  similar  purpose. 


III.  NEW  CONSTITUTIVE  MODELING 
APPROACHES  EOR  NEAR-|0  0  1|  CREEP 
IN  THE  LH  REGIME 

Before  describing  in  detail  the  new  constitutive  model, 
it  is  important  to  clarify  the  representative  microstruc- 
tural  features  on  which  the  new  creep  model  is  based. 
Various  microstructural  and  defect  features  at  different 
length  scales  may  limit  thermomechanical  responses  of 
single-crystal  superalloys.  Figure  3  shows  typical  micro¬ 
structures  of  a  single-crystal  superalloy  viewed  at  differ¬ 
ent  length  scales  (magnifications).  At  the  macroscopic 
scale  (several  millimeters  to  centimeters,  Figure  3(a)) 
columnar  dendrites,  which  are  slightly  misoriented 
relative  to  each  other,  and  a  colony  of  a  few  misoriented 
grains  (fleckles)  may  be  influential  features.  At  the 
intermediate  scale  (about  a  few  hundred  micrometers. 
Figure  3(b)),  the  chemistry  gradients  between  dendritic 
cores  and  the  interdendritic  region  that  contains  non- 
metallic  inclusions  and  microvoids  are  probably  key 
features  that  one  needs  to  account  for.  At  the  micro¬ 
scopic  scale  (several  to  tens  of  micrometers.  Figure  3(c)), 
regularly  distributed  cuboidal  7'  precipitates  and  sur¬ 
rounding  7-matrix  channels  are  microstuctural  features 
that  control  thermomechanical  responses.  The  new  creep 
constitutive  model  was  developed  by  incorporating 
dislocation  micromechanisms  identified  from  electron 
microscopy  analysis^^  together  with  certain  assump¬ 
tions  about  strain  hardening.  In  this  sense.  Figure  3(c) 
represents  the  microstructural  features  that  are  repre¬ 
sented  by  the  current  creep  constitutive  model.  The  new 
model  has  no  explicit  treatment  of  the  other  microstruc¬ 
tural  and  defect  features  shown  in  Figures  3(a)  and  (b). 
However,  dislocation  behaviors  in  the  7/7'  microstruc¬ 
ture  (Figure  3(c))  seem  to  be  a  major  factor  that  controls 
anisotropic  primary  and  secondary  creep  in  the  LH 
regime,  which  is  the  focus  of  the  current  study. 

In  the  new  constitutive  model,  the  plastic  shear  strain 
rate  caused  by  creep  deformation  consists  of  two 
contributions 

r=fyyl  +  yl^  [1] 

where  7|j,  is  the  shear  strain  rate  caused  by  the  a/2(110) 
dislocation  activity  (the  slip  system  index  a  ranges  from  1 
to  12)  in  7-matrix  channels  and/y  is  the  volume  fraction  of 
the  7  matrix.  Also,  7gp  is  the  shear  strain  rate  caused  by 
shearing  7'  precipitates  and  the  7  matrix  by  a(112) 
dislocation  ribbons  (a  ranges  from  13  to  24).  Details 
about  how  constitutive  formulations  were  developed  for 
each  of  the  two  contributions  are  given  in  the  subsequent 
sections. 

A.  7^  Resulting  from  a/2(110)  Dislocations  in  the  7 
Matrix 

As  noted  previously,  7“  reflects  a/2(110)  dislocation 
behaviors  in  narrow  7-matrix  channels.  The  Orowan 
equation  was  used  to  express  7“ 


Fig.  3  Microstructural  and  defect  features  that  can  influence  ther 
momechanical  behaviors  of  single  crystal  superalloys  at  (a)  macro 
scopic,  (/;)  intermediate,  and  (c)  microscopic  length  scales.  All 
micrographs  were  taken  from  a  second  generation,  single  crystal 
superalloy  PWA1484.  The  directional  solidification  is  along  the  verti 
cal  direction  in  (a),  whereas  (b)  and  (c)  are  the  cross  sectional  view. 


yl  =  bplvl  [2] 

where  b,  p“,  and  v“  are  the  magnitude  of  Burgers 
vector,  the  density  of  a/2(110)  dislocations  in  the  7 
matrix,  and  the  corresponding  average  velocity,  respec¬ 
tively.  In  Eq.  [2],  the  evolution  of  was  assumed  to 
follow  the  Kocks-Mecking-Estrin  equation  (strain 
hardening  law)^^^’^"*'  by 
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Pm  =  (ki\/p“  -k2P^)|)'^|  [3] 

where  ki  and  k2  are  the  constants  that  control  the  con¬ 
tributions  of  dislocation  multiplication  (hardening) 
and  annihilation  (recovery),  respectively,  for  the  evolv¬ 
ing  a/2(110)  dislocation  density.  In  Eq.  [2],  a  power 
law  was  adopted  to  describe 

a  1 

<  =  sign(4)  [4] 

sm 

and 

».  =  <exp(-|;)  |5| 

where  and  g„,  are  the  resolved  shear  stress  and  the 
slip  resistance  for  a/2(110)  dislocations  in  the  y  matrix, 
respectively,  and  m  is  the  strain-rate  sensitivity  param¬ 
eter.  Also,  in  Eq.  [5],  and  Q,„  are  the  preexponent 
and  a  phenomenological  activation  energy  for  time- 
dependent  a/2(110)  dislocation  motion  in  the  y  matrix, 
respectively.  In  this  equation,  the  slip  resistance  g„,  for 
a/2(110)  dislocations  in  the  y  matrix  was  set  to  consist 
of  three  contributions 


Sm  —  So  “f  .?Orowaii  “f  ^/j7SF  [h] 

where  go,  gorowan,  and  g^sF  are  the  slip  resistances 
caused  by  the  intrinsic  y-matrix  strength,  Orowan  bow¬ 
ing  in  the  y-matrix  channels,  and  the  interaction  of 
a/2(110)  dislocations  with  a(112)  dislocation  ribbons, 
respectively.  The  Orowan  bowing  contribution  gorowan 
was  expressed  by  rj'fih/l,  where  t]',  and  /  are  the  geo¬ 
metric  constant,  the  shear  modulus,  and  the  width  of 
the  y-matrix  channel  in  the  {111}  plane  (/  =  (3/2) 
where  ly  is  the  channel  thickness),  respectively.  In 
Eq.  [6],  the  evolution  of  g,„sF  was  expressed  by 

Lsf  =  Iir=i3  ''-sfTsf  [7] 

and 

(i+i:/l„Asp(|;))  IS] 

where  h,,„  and  are  the  initial  hardening  modulus 
and  the  coefficient  for  the  a/2(110)  -  a(112)  interaction, 
respectively.  Also,  y^p,  t^^p,  and  g^^p  are  the  shear  rate, 
the  resolved  shear  stress  and  the  slip  resistance  for  the 
(112){1  1  1}  slip  system  respectively.  Details  about 
these  three  state  variables  will  be  discussed  in  the  follow¬ 
ing  section.  Equations  [7]  and  [8]  were  introduced  to 
account  for  strengthening  of  the  a/2(110)  slip  systems 
because  of  the  interaction  with  a(112)  dislocation 
ribbons. 


B.  The  a(112)  Dislocation  Shear  Rale,  yjp 

Similar  to  Eq.  [2],  ygp  was  expressed  using  the 
Orowan  equation 

fSF  =  ^PSF^'SF  P] 


where  p^p  and  v^p  are  the  density  of  a(112)  dislocation 
ribbons  and  the  corresponding  average  velocity, 
respectively.  In  Eq.  [9],  the  pgp  evolution  (p^p)  and  vgp 
were  described  numerically  by  adopting  dislocation 
evolution  equations  for  the  nucleation,  propagation, 
and  hardening  of  a(112)  SF  dislocation  ribbons,  as 
proposed  by  Rae  et  The  evolution  of  p|p  was 

expressed  by 


[10] 


where  k^Po,  ksp],  and  ksF2  are  the  constants  that  control 
the  contribution  of  the  corresponding  tenns.  Equation  [10] 
basically  incorporates  the  nucleation,  propagation,  and 
immobilization  of  a(112)  dislocation  ribbons. 

The  first  term  of  the  right-hand  side  (RHS)  of  Eq.  [10] 
accounts  for  the  production  of  the  a(112)  dislocation 
density  because  of  the  nucleation  of  a(112)  dislocation 
ribbons  from  the  reaction  of  two  a/2(110)  disloca¬ 
tions  in  the  y  matrix  and  the  y/y'  interface. Table  I 
shows  a  list  of  slip  systems  for  a/2(110)  slip  and  a(112) 
slip  along  with  the  corresponding  notations. 

The  reaction  of  two  a/2(110){l  1  1}  slip  systems  can 
produce  the  a  (112){1  1  1}  slip  system  if  the  resulting 
(112)  vector  and  (110)  vector  lie  in  the  same  {111} 
plane.  Based  on  this  reaction  criterion.  Table  II  lists 
possible  a(l  12){1  1  1}  slip  systems  from  reactions  of  two 
a/2(110){l  1  1}  slip  systems. 

Because  the  reaction  matrix  in  Table  II  is  symmetric, 
only  a  diagonal  half  of  the  reaction  matrix  needs  to  be 
considered.  Elere,  one  can  notice  from  Table  II  that  each 
a(112){l  1  1}  slip  system  can  be  produced  by  three 
different  pairs  of  a/2(110){l  1  1}  slip  systems.  This 
finding  implies  that  the  availability  of  nucleating  each  of 
a(112){l  1  1}  slip  system  can  be  checked  by  examining 
activities  of  the  corresponding  three  pairs  of  a/2(110){l 
1  1}  slip  systems.  Thus,  the  nucleation  term  in  Eq.  [10]  is 
completed  by  introducing  Mgp  to  quantify  the  avail¬ 
ability  of  the  a(112){l  1  1}  nucleii 

=  max 

\  i=  1,2,3 

[11] 


max(a,-,6,-) 


avg(ai,*i)  I  I  _1 _ jyj  ] 


Table  I.  Twelve  a/2(110){l  1  1}  and  Twelve  a(112){l  1  1} 
Slip  Systems  and  Their  Notations 


(110){1 

1 1} 

(112){1  1  1} 

A1 

oil 

(111) 

A  I  ]112](111) 

A2 

101 

(111) 

A  II  ]121](111) 

A3 

110 

(111) 

A  III  ]211](111) 

B1 

Oil 

(111) 

BI  [112]  (111) 

B2 

101 

111 

B  II  [121]  (111) 

B3 

110] 

111) 

Bill  [2111(111) 

Cl 

Oil] 

111) 

Cl  [112]  (111 

C2 

101 

111 

C  II  [121](111) 

C3 

110 

111 

C  III  [211]  (111) 

D1 

[Oil 

111 

D  I  [112](111) 

D2 

lOf 

(111) 

D  II  [1211(111) 

D3 

]110](111) 

D  III  [211]  (111) 
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Table  II.  The  Matrix  Showing  Output  a(112){l  1  1}  Slip  Systems  from  Reactions  of  Two  a/2(110){l  1  1}  Slip  Systems 


A1 

A2 

A3 

B1 

B2 

B3 

Cl 

C2 

C3 

D1 

D2 

D3 

A1 

A  I 

A  II 

B  I 

B  II 

A  I 

A  II 

A2 

A  I 

A  III 

A  I 

C  I 

C  III 

A  III 

A3 

A  II 

A  III 

A  II 

A  III 

D  II 

D  III 

B1 

A  I 

A  II 

B  I 

B  II 

B  II 

B  I 

B2 

B  I 

B  I 

B  III 

B  III 

D  I 

D  III 

B3 

B  II 

B  II 

B  III 

C  II 

C  III 

B  III 

Cl 

C  I 

C  II 

C  I 

C  II 

D  I 

D  II 

C2 

A  I 

A  III 

C  III 

C  I 

C  III 

C  I 

C3 

C  III 

B  II 

B  III 

C  II 

C  III 

C  II 

D1 

D  II 

D  I 

C  I 

C  II 

D  I 

D  II 

D2 

D  III 

B  I 

B  III 

D  I 

D  I 

D  III 

D3 

A  II 

A  III 

D  III 

D  II 

D  II 

D  III 

where  for  each  of  a(112){l  1  1}  slip  systems  (a  =  13  to 
24  in  Eq.  [11]),  the  index  i  indicates  three  available 
pairs  of  a/2(110){l  1  1}  slip  systems  and 

and  are  the  smaller,  larger,  and 

average  resolved  shear  stresses  for  each  a/2(110){l  1  1} 
pair,  respectively.  Also,  tsf  is  the  stress  required  to 
create  a  stacking  fault  (tsf  =  >7sF-^SF/h,  where  rjsF  and 
Esf  are  the  constant  and  the  stacking  fault  energy  in 
the  7  matrix,  respectively),  and  and  are 

constants.  The  first  term  of  the  RHS  of  Eq.  [11]  was 
formulated  such  that  the  nucleation  of  the  a(112){l  1 
1}  comes  from  a  most  active  a/2(110){l  1  1}  pair  by 
checking  the  corresponding  interactions  in  Table  II.  In 
Eq.  [II],  the  second  term  Mq  was  intended  to  cut  off 
the  value  range  of  the  first  term  before  the  normaliza¬ 
tion  by  R,v/. 

The  second  term  of  the  RHS  of  Eq.  [10]  is  intended  to 
incorporate  the  multiplication  of  a(112)  dislocation 
ribbons.  Here,  /sf  is  the  average  y'  SF-shearing  distance 
(/sF  =  where  L  is  the  size  of  the  cuboidal 

y'  precipitate^^^').  The  third  term  of  the  RHS  of  Eq.  [10] 
accounts  for  the  immobilization  of  a(112)  dislocation 
ribbons  by  interactions  between  different  a(112)  dislo¬ 
cation  ribbons.  The  law  governing  the  immobilization 
follows  a  modified  Kock,  Argon,  and  Ashby  form  that 
accounts  for  the  maximally  stressed  a(112){l  1  1}  glide 
system.  Here,  T/^f  expressed  by 


=  1  —  exp 


(max 

/?=13...24 

/J/a 


1  fP 
p/SE 

^SF 

-max 

1 

^SF 

1 

I  .^niax  I 
I  ^SF  I 


+  PSF 


qSF' 


avg(a,j8) 

‘'SF 


[12] 


In  Eq.  [12],y^p  is  the  interaction  coefficient,  and  psr 
and  qsF  are  constants.  Also,  1™“  is  the  maximum 
resolved  shear  stress  of  all  a(112){l  1  1}  slip  systems, 
and  is  the  average  resolved  shear  stress  of  the 

two  a(112){l  1  1}  slip  systems.  Equation  [12]  was 
formulated  in  such  a  fashion  so  that  rapid  immobiliza¬ 
tion  (in  a  form  of  an  S-type  curve)  of  the  primary  slip 
system  takes  place  as  the  secondary  slip  system  becomes 
active. 


Returning  to  Eq.  [9],  the  average  velocity  of  a(112) 
dislocation  ribbons  was  expressed  using  a  power  law  as 


VsF  =  '"SFo 


■•SF 


^SF 


1/msF 


sign(T“F) 


[13] 


and 


i'SFo  =  v'sp„exp(^-^^  [14] 

where  nrsF  is  the  strain  rate  sensitivity  parameter. 
Also,  in  Eq.  [14],  v'sp„  and  Qsf  are  the  preexponent 
and  the  activation  energy  for  a(112)  dislocations, 
respectively.  Here,  the  slip  resistance  g^p  for  a(112) 
dislocation  ribbons  was  set  to  originate  from  two  con¬ 
tributions 


The  first  and  second  terms  of  the  RHS  of  Eq.  [15]  are 
hardening  contributions  from  a(112)  -  a(112)  interac¬ 
tions  and  from  a(112)  -  a/2(110)  interactions,  respec¬ 
tively.  The  initial  value  of  g^p  was  determined  by 
isF  =  -^APB /h,  where  Eapb  is  the  APB  energy  for  y' 
precipitates.  In  Eq.  [15],  7™“  is  the  maximum  shear 
rate  of  a(112){l  1  1}  slip  systems.  Here,  /i^f  is  the 
hardening  modulus  from  a(112)  -  a(112)  interactions 
and  was  expressed  by 

hip  =  ■  sech'/2  ^  •  (1  +  ksf^^f)  [16] 

where  /?,  and  7;  are  the  initial  hardening  modulus  and 
the  normalization  constant,  respectively.  Also,  «sf  is 
the  scaling  constant  that  links  Hip  of  Eq.  [12]  to  hip. 
For  the  second  term  in  Eq.  [15],  /7“f  is  the  hardening 
modulus  stemming  from  a(112)  -  a/2(110)  interactions 
that  was  expressed  by 

=  /^SF.  •  (1  +  tanh  J,)  )  [17] 

where  h^p,„  and  are  the  initial  hardening  modulus 
and  the  coefficient  for  a(112)  -  a/2(110)  interactions, 
respectively. 
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IV.  SIMULATION  OF  ANISOTROPIC  CREEP 
FOR  CMSX-4 

The  new  constitutive  model  described  in  the  previous 
section  was  used  to  capture  orientation  sensitivity  of 
near-[0  0  1]  creep  of  a  single-crystal  superalloy  CMSX-4. 
All  constitutive  formulations  (Eqs.  [1]  to  [17])  were 
implemented  in  the  commercial  FEM  package  ABA- 
QUS  (Simulia,  Providence,  RI)  through  the  user  mate¬ 
rial  subroutine  (UMAT).  Here,  a  tangent  modulus 
methodP^’l  was  adopted  for  the  time  integration. 
Preliminary  parametric  studies  were  performed  to  iden¬ 
tify  reasonable  ranges  for  values  of  major  input  param¬ 
eters.  Table  III  lists  values  of  input  parameters  for  the 
description  of  7“,  and  y^p. 

In  Table  III,  the  values  of  parameters  marked  with 
the  asterisk  were  determined  by  fitting  to  the  creep  data 
in  Figure  5.^^^  The  values  for  activation  energies  (Q„,  and 
gsp)  and  fault  energies  (Eapb  and  Esf)  were  approxi¬ 
mated  from  the  experimental  and  modeling  data.^^^ 
The  values  for  ly,  fy,  L  and  b  were  directly  taken  from 
the  input  microstructure  (cuboidal  y'  precipitates  with 
an  edge  length  of  0.52  {.iva  and  the  volume  fraction  of 
68  pet).  The  constants  in  Eqs.  [3]  and  [10],  and  the 
parameters  in  Eqs.  [4]  and  [13]  were  determined  from  a 
range  of  values  used  in  other  models. The  geomet¬ 
ric  constants  r]'  and  were  set  as  1  for  the  current 
simulations.  Also,  the  coefficients  and  /sf)  for 

the  interaction  of  dilferent  slip  mechanisms  were 
assumed  to  be  1  because  their  values  are  not  informed. 
The  shear  modulus  /r  was  calculated  from  /r  = 
[0.5'C44'(Cii  -  Ci2)]’^^,  where  the  temperature  depen¬ 
dence  of  elastic  constants  (Cn,  C12,  and  C44)  for  CMSX- 
4  was  determined  by  experimental  fits  from  Allan. 
Figure  4(a)  shows  the  sample  geometry  used  for  FEM 
simulations.  A  cylindrical  sample  was  meshed  by  the 
combination  of  linear  wedge  (C3D6)  and  brick  (C3D8) 
elements.  Boundary  conditions  for  creep  deformation 
were  applied  to  the  meshed  sample  by  restraining  the 
bottom  face  vertically  and  applying  appropriate  body 
forces  to  the  top  mesh  layer  of  the  cylindrical  sample  to 
induce  an  axial  stress  comparable  with  the  creep  stress. 
Figure  4(b)  shows  the  contour  of  accumulated  shear  in 
the  deformed  sample  geometry  after  15  pet  creep  strain 
along  the  single  slip-oriented  direction.  One  can  observe 


Table  III.  Values  of  Input  Parameters  Used  for  the  Current 
Simulations 


ACC 

1  m 

fsp 

Qm 

323  kJ/mol 

2sf 

280  kJ/mol 

M„ 

0.4* 

go 

300  MPa* 

Aapb 

0.12  J/nf 

Rm 

0.3* 

ki 

1.2  X  10*/m 

Esf 

0.11  J/m^ 

h. 

40  GPa* 

kl 

30 

ksFo 

4 

IT 

0.07* 

< 

1  X  10®/s 

ksFl 

2  X  lO^"* 

WSF 

24* 

m 

0.13 

ksF2 

2 

hsFm 

10  GPa* 

n' 

1 

r’sFo 

1  X  lO^/s 

1 

ly 

0.068  X  10^®  m 

msF 

0.25 

/SF 

1 

h 

'^sm 

10  GPa* 

bSF 

1 

PSF 

0.028* 

/mSF 

fy 

b 

1 

0.32 

2.489  A 

L 

ksFo 

0.52  X  10^®  m 
4 

?SF 

250* 

that  the  deformed  geometry  is  locally  nonuniform, 
which  is  typical  for  single  crystals  deformed  in  low 
symmetry  {viz.,  single  slip  oriented)  orientations. 

Using  the  new  constitutive  model  described  in  Section 
III  and  the  sample  geometry  shown  in  Figure  4(a),  the 
creep  simulations  were  performed  in  three  dilferent 
crystallographic  directions  near  the  [0  0  1]  at  1023  K 
(750  °C)  and  750  MPa,  and  the  resulting  creep  curves 
were  compared  with  those  from  CMSX-4^^^^  in  Figure  5. 
Here,  three  creep  orientations  K,  L,  and  M  are 
misoriented  from  the  [0  0  1]  by  5.3,  8.0,  and  17.5  deg, 
respectively.  The  prediction  result  (Figure  5)  confirms 
that  the  new  creep  constitutive  model  reasonably  the 
orientation  sensitivity  of  primary  and  secondary  creep 
captures  for  near-[0  0  1]  oriented  single-crystal  superal¬ 
loys.  In  Figure  6,  the  shear  strain  rates  from  a(112)  slip 
(y^p)  and  a/2(110)  matrix  slip  (}>,„)  were  averaged  over 
an  entire  sample  volume  and  plotted  as  a  function  of 
creep  time.  As  expected,  a(112)  slip  dominates  for  all 
three  orientations,  and  the  degree  of  primary  creep  is 
shown  to  be  inversely  proportional  to  the  level  of  the  y^p 
suppression,  which  is  highest  in  M  and  lowest  in  K. 
Here,  it  is  interesting  to  note  that  y,„  shows  the  opposite 
trend,  the  highest  suppression  in  K  and  the  lowest 
suppression  in  M.  Importantly,  this  result  implies  that 
a/2(110)  matrix  slip  cannot  be  ignored,  and  it  plays  an 
important  role  in  blocking  a(112)  dislocation  ribbons 
slip  activities  effectively,  leading  to  the  transition  to 
secondary  creep.  The  current  simulation  results  seem  to 
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Fig.  4  (a)  A  cylindrical  sample  geometry  meshed  by  linear  wedge 

(C3D6)  and  brick  (C3D8)  elements,  (h)  The  contour  of  accumulated 
shear  in  the  deformed  sample  geometry  after  15  pet  creep  strain  in  a 
single  slip  oriented  direction. 
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Fig.  5  Creep  curve  comparison  between  predictions  (solid  curves) 
and  CMSX  4  data  (circles):  three  creep  loading  directions  K,  L,  and 
M  are  indicated  in  a  crystallographic  triangle.  Experimental  creep 
curves  were  taken  from  Ref.  3. 
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0  20  40  60  80  100 


Fig.  6  Average  shear  strain  rates  from  a(112)  SF  slip  (ygp,  solid 
lines)  and  a/2(110)  matrix  slip  (y^^,  broken  lines)  plotted  as  a  func 
tion  of  creep  time  for  orientations  K,  L,  and  M. 

be  in  agreement  with  conclusions  of  experimental  studies 
that  emphasized  the  importance  of  a/2(110)  matrix 
dislocations  to  control  a(112)  slip  propagation."*’® 


V.  CONCLUDING  REMARKS 

The  current  study  developed  a  new  crystallographic 
constitutive  model  to  predict  highly  anisotropic  primary 
and  secondary  creep  near  the  [0  0  1]  orientation  of 
single-crystal  superalloys  in  relatively  low-temperature 
and  high-stress  ranges.  Crystal  plasticity  formulations 
using  slip  systems  kinematics  were  developed,  guided  by 
experimentally  observed  dislocation  micromechanisms. 
The  new  constitutive  model  captured  the  experimental 
creep  behavior  of  near-[0  0  1]  oriented  CMSX-4 


reasonably.  It  suggests  that  crystal  plasticity  is  a 
reasonable  framework  that  can  be  used  to  capture 
primary  and  secondary  creep  of  single-crystal  superal¬ 
loys  in  the  low-temperature  and  high-stress  regime.  The 
simulation  results  also  suggest  that  in  addition  to  a(112) 
-  a(112)  interactions,  a/2(110)  matrix  dislocations  are 
responsible  for  limiting  the  propagation  of  a(112) 
dislocation  ribbons  by  a(112)  -  a/2(110)  interactions. 
This  was  in  agreement  with  the  suggestion  made  by 
experimental  studies. 
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